Two electrons on a hypersphere: a quasi-exactly solvable model 
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We show that the exact wave function for two electrons, interacting through a Coulomb potential 
but constrained to remain on the surface of a D-sphere (T> > 1), is a polynomial in the interelectronic 
distance u for a countably infinite set of values of the radius R. A selection of these radii, and the 
associated energies, are reported for ground and excited states on the singlet and triplet manifolds. 
We conclude that the T> — 3 model bears the greatest similarity to normal physical systems. 
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Quantum mechanical models for which it is possible to 
solve explicitly for a finite portion of the energy spectrum 
are said to be quasi-exactly solvable [l[ . They have ongo- 
ing value and are useful both for illuminating more com- 
plicated systems and for testing and developing theoreti- 
cal approaches, such as density functional theory (DFT) 
and explicitly correlated methods [BHU . One of the 
most famous two-body models is the Hooke's law atom 
which consists of a pair of electrons, repelling Coulombi- 
cally but trapped in a harmonic external potential with 
force constant k. This system was first considered nearly 
50 years ago by Kestner and Sinanoglu H, s olved ana- 
lytically in 1989 for one particular k value [lOj, and later 
for a countably infinite set of k values (llj . 

A related system consists of two electrons trapped on 
the surface of a sphere of radius R. This has been used 
by Berry and collaborators [l2-15| to understand both 



weakly and strongly correlated systems and to suggest an 
"alternating" version of Hund's rule [l6|. Seidl utilized 
this system to develop new correlation functionals (l7j 



within the adiabatic connection in DFT [18|. We will 
use the term "sphcrium" to describe this system. 

In recent work [l9j . we examined various schemes and 
described a method for obtaining near-exact estimates of 
the 1 S ground state energy of spherium for any given R. 
Because the corresponding Hartree-Fock (HF) energies 
are also known exactly [19J , this is now one of the most 
complete theoretical models for understanding electron 
correlation effects. 

In this Letter, we consider 2?-spherium, the generaliza- 
tion in which the two electrons are trapped on a 2?-sphcrc 
of radius R. We adopt the convention that a P-sphcre 
is the surface of a {T> + l)-dimensional ball. (Thus, for 
example, the Berry system is 2-spherium.) We show that 
the Schrodinger equation for the 1 S and the 3 P states can 
be solved exactly for a countably infinite set of R values 
and that the resulting wave functions are polynomials in 
the interelectronic distance u = \ri — r%\- Other spin and 
angular momentum states can be addressed in the same 
way using the ansatz derived by Breit [20j | . 



The electronic Hamiltonian, in atomic units, is 
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and because each electron moves on a T>- sphere, it is 
natural to adopt hyperspherical coordinates [2l|, |22J . 

For 1 S states, it can be then shown [l9[ that the wave 
function S(u) satisfies the Schrodinger equation 
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By introducing the dimensionless variable x = u/2R, this 
becomes a Heun equation [23( with singul ar p oints at 
x = —1, 0, +1. Based on our previous work [19( and the 
known solutions of the Heun equation [24j , we seek wave 
functions of the form 



S{u) = ^ S k ' 



(3) 



k=0 



and substitution into <j2j) yields the recurrence relation 
Sfc+i + [Hk + 2V-2) 7 ^~E]s k 
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(k + 2)(k + V) 



with the starting values 



{so,«i} = 



{0,1} V = l 

{l,l/(2?-l)} V>2 

Thus, the Kato cusp conditions [25| are 

for electrons on a circle (T> = 1) and 
S'(0) 1 
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in higher dimensions. We note that the "normal" Kato 
value of 1 /2 arises for T> = 3, suggesting that this may the 
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TABLE I: Radius R, energy E and wave function S(u) or T(u) 
of the first 1 S and 3 P polynomial solutions for two electrons 
on a D-sphere 



State 


V 


27? 


E 


S(u) or T(u) 




1 


Ve 


2/3 


u(l + w/2) 


KS 


2 




1 


1 + u 




3 


Vw 


1/2 


l + u/2 




4 




1/3 


l + w/3 




1 


V6 


1/2 


l+u/2 


3p 


2 


\/l5 


1/3 


l + u/3 




3 


728 


1/4 


l + w/4 




4 


%/45 


1/5 


1 + it/5 



most appropriate model for atomic or molecular systems. 
We will return to this point below. 

The wave function (J3J) reduces to the polynomial 



S n ,m (U) 



(8) 



(where m the number of roots between and 2R) if, 
and only if, s„ + i — s n+ 2 = 0. Thus, the energy E n . m 
is a root of the polynomial equation s n +i = (where 
degSn+i = l(n + 1)/2J) and the corresponding radius 
Rn,m is found from ([4]) which yields 
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S n .m(u) is the exact wave function of the m-th excited 
state of S symmetry for the radius Rn, m - 
If we write the 3 P state wave function as [2(1 



? * = (cos#i - cos 6 2 )T(u) 



(10) 



where 9\ and 9 2 are the 2?-th hyperspherical angles of 
the two electrons [U 22 1, the symmetric part satisfies 
the Schrddinger equation 
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and the antisymmetric part provides an additional ki- 
netic energy contribution V/(2R 2 ). 

Substituting the power series expansion 



T{u) 



fc=0 



into (fTT|) yields the recurrence relation 

t k +i + [k(k + W)-jp-E] t k 



tk+2 — 



V + 2) 



[k + 2){k 
with the starting values 

{t ,ti} = {l,l/(Z> + l)} 
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yielding the cusp condition 
T'(0) 



T(0) V+l 
The wave function (|T2")) reduces to the polynomial 
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when the energy E n<m is a root of t n+ \ = and the cor- 
responding radius R n ,m is found from (|13p which yields 
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T n .m(u) is the exact wave function of the m-th excited 
state of 3 P symmetry for the radius R niTn - 

It is illuminating to begin by examining the simplest 
1 5 and 3 P polynomial solutions. Except in the V = 1 
case, the first *S solution has 



/ (2P-l)(2P-2) 

-n-1,0 = \/ ^ -^1.0 
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and the first 3 P solution has 



(2X>+l)(22? + 2) 



£■1,0 — 



V+l 
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(19) 



These are tabulated for V = 1,2,3,4, together with the 
associated wave functions, in Table [IJ 

In the V = 1 case (i.e. two electrons on a circle), the 
first singlet and triplet solutions have £2,0 = 2/3 and 
£4,0 = 1/2, respectively, for the same value of the radius 
(Vd/2 w 1.2247). The corresponding wave functions are 
related by 52,0 = uTi^. Unlike T1.0, the singlet wave- 
function 52,0 vanishes at u = 0, and exhibits a second- 
order cusp condition, as shown in ([6]). 

For the 2-spherium (V = 2 case), we know from our 
previous work [l9[ that the HF energy of the lowest 1 S 
state is T^hf — 1/R. It follows that the exact correlation 
energy for R = V3/2 is E COII = 1 - 2/V3 -0.1547 
which is much larger than the limiting correlation ener- 
gies of the helium-like ions (—0.0467) [26( or Hooke's law 
atoms (—0.0497) [27j. This confirms our view that elec- 
tron correlation on the surface of a sphere is qualitatively 
different from that in three-dimensional physical space. 

The 3-spherium (V = 3 case), in contrast, possesses 
the same singlet and triplet cusp conditions — Eqs. 
(|7|) and (fl~5|) — as those for electrons moving in three- 
dimensional physical space. Indeed, the wave functions 
in Tabled 



Si,o(«) = l + u/2 
T h0 (u) = l + u/4 



(R=V5/Z) 
(R = V7) 



(20) 
(21) 



have precisely the form of the ansatz used in Kutzclnigg's 
increasingly popular R12 methods 0, Q. Moreover, it 
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TABLE II: Radii R n>m and energies E n , m for 1 S states of two electrons on a D-sphere (25=1,2,3) 









V 


= 1 






V 


= 2 






V 


= 3 






n/m 





1 


2 


3 





1 


2 


3 





1 


2 


3 




1 










0.8660 








1.5811 










2 


1.2247 








2.6458 








4.0620 










3 


3.3912 








5.4312 


1.4150 






7.5154 


2.2404 








4 


6.5439 


1.9178 






9.2211 


3.7379 






11.961 


5.3320 






- 

K 


5 


10.693 


4.7071 






14.012 


7.0848 


1.9256 




17.404 


9.3775 


2.8554 






6 


15.841 


8.4583 


2.5522 




19.804 


11.448 


4.7683 




23.846 


14.410 


6.5350 






7 


21.989 


13.199 


5.9404 




26.597 


16.817 


8.6593 


2.4123 


31.287 


20.439 


11.158 


3.4415 




8 


29.136 


18.936 


10.277 


3.1515 


34.389 


23.190 


13.583 


5.7566 


39.728 


27.466 


16.768 


7.6903 




1 










1.0000 








0.5000 










2 


0.6667 








0.2857 








0.1818 










3 


0.1957 








0.1271 


1.8729 






0.0930 


1.0459 






a 


4 


0.0934 


1.0875 






0.0706 


0.4294 






0.0559 


0.2814 






a 
w 


5 


0.0547 


0.2821 






0.0446 


0.1743 


2.3597 




0.0371 


0.1279 


1.3798 






6 


0.0359 


0.1258 


1.3817 




0.0306 


0.0916 


0.5278 




0.0264 


0.0722 


0.3512 






7 


0.0253 


0.0703 


0.3471 




0.0223 


0.0557 


0.2100 


2.7065 


0.0197 


0.0461 


0.1546 


1.6253 




8 


0.0188 


0.0446 


0.1515 


1.6110 


0.0169 


0.0372 


0.1084 


0.6035 


0.0152 


0.0318 


0.0854 


0.4058 



TABLE III: Radii R n ,m and energies E n>m for 3 P states of two electrons on a D-sphere (D=l,2,3) 
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= 1 






V 


= 2 






V 


= 3 






n/m 





1 


2 


3 





1 


2 


3 





1 


2 


3 




1 


1.2247 








1.9365 








2.6458 










2 


3.3912 








4.7958 








6.2048 








CO 


3 


6.5439 


1.9178 






8.6227 


2.6738 






10.718 


3.4111 






■3 


4 


10.693 


4.7071 






13.435 


6.2041 






16.205 


7.6748 






oi 


5 


15.841 


8.4583 


2.5522 




19.241 


10.665 


3.3588 




22.678 


12.852 


4.1285 






6 


21.989 


13.199 


5.9404 




26.043 


16.094 


7.5340 




30.142 


18.979 


9.0701 






7 


29.136 


18.936 


10.277 


3.1515 


33.842 


22.505 


12.615 


4.0095 


38.600 


26.077 


14.897 


4.8130 




8 


37.283 


25.671 


15.599 


7.1177 


42.640 


29.907 


18.650 


8.8083 


48.054 


34.155 


21.654 


10.411 




1 


0.5000 








0.3333 








0.2500 










2 


0.1739 








0.1304 








0.1039 










3 


0.0876 


1.0196 






0.0706 


0.7343 






0.0588 


0.5801 






SB 

o 


4 


0.0525 


0.2708 






0.0443 


0.2078 






0.0381 


0.1698 






H 


5 


0.0349 


0.1223 


1.3433 




0.0304 


0.0989 


0.9972 




0.0267 


0.0832 


0.8067 






6 


0.0248 


0.0689 


0.3401 




0.0221 


0.0579 


0.2643 




0.0198 


0.0500 


0.2188 






7 


0.0186 


0.0439 


0.1491 


1.5858 


0.0168 


0.0380 


0.1210 


1.1974 


0.0153 


0.0335 


0.1025 


0.9821 




8 


0.0144 


0.0303 


0.0822 


0.3948 


0.0132 


0.0268 


0.0690 


0.3093 


0.0121 


0.0240 


0.0597 


0.2583 



can be shown [28[ that, as R — > 0, the correlation energy 
Ecoti approaches —0.0476, which nestles between the cor- 
responding values for the helium-like ions (—0.0467) 
and the Hooke's law atom (-0.0497) [13 . Again, this 
suggests that the T> = 3 model ("electrons on a hyper- 
sphere") bears more similarity to common physical sys- 
tems than the T> = 2 model ( "electrons on a sphere" ) . 

Numerical values of the energies and radii, for polyno- 
mial wave functions in T> = 1, 2, 3, are reported in Table 



|II](for l S states) and Table [Hi] (for 3 P states). 

For fixed T>, the radii increase with n but decrease 
with m, and the energies behave in exactly the opposite 
way. As R (or equivalently n) increases, the electrons 
tend to localize on opposite sides of the sphere, a phe- 
nomenon known as Wigner crystallization 29| which has 
also been observed in other systems [H|,|30|. As a result, 
for large R, the ground state energies of both the singlet 
and triplet state approach 1/(2R). Analogous behavior 
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is observed when V — > oo 31, 3 



In conclusion, we have shown that the system of two 
electrons, interacting via a Coulomb potential but con- 
strained to remain on a D-sphere, can be solved exactly 
for an infinite set of values of the radius R. We find 
that the 3-spherium (D = 3 model), wherein the elec- 
trons are confined to a three-dimensional surface of a 
four-dimensional ball, has greater similarity to normal 
physical systems than the more familiar D = 2 case. 

We believe that our results will be useful in the future 
development of correlation functionals within density- 



functional theory [33[, intracule functional theory [3J- 
IHsl, and explicitly correlated methods [B-Q. They also 
shed new light on dimension-dependent correlation ef- 
fects, and may be used as an alternative system for study- 
ing quantum dots |40j |. 
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